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ABSTRACT 

Cold massive cores are one of the earliest manifestations of high mass star formation. 
Following the detection of SiO emission from G333. 125-0.562, a cold massive core, 
further investigations of the physics, chemistry and dynamics of this object has been 
carried out. Mopra and NANTEN2 molecular line profile observations, Australia Tele- 
scope Compact Array (ATCA) line and continuum emission maps, and Spitzer 24 and 
70 /im images were obtained. These new data further constrain the properties of this 
prime example of the very early stages of high mass star formation. A model for the 
source was constructed and compared directly with the molecular line data using a 
3D molecular line transfer code - MOLLIE. The ATCA data reveal that G333.125- 
0.562 is composed of two sources. One of the sources is responsible for the previously 
detected molecular outflow and is detected in the Spitzer 24 and 70 /im band data. 
Turbulent velocity widths are lower than other more active regions of G333 which 
reflects the younger evolutionary stage and/or lower mass of this core. The molecu- 
lar line modelling requires abundances of the CO isotopes that strongly imply heavy 
depletion due to freeze-out of this species onto dust grains. The principal cloud is 
cold, moderately turbulent and possesses an outflow which indicates the presence of a 
central driving source. The secondary source could be an even less evolved object as 
no apparent associations with continuum emissions at (far-)infrared wavelengths. 
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1 INTRODUCTION 

High mass star formation is unlikely to be simply a scaled 
up version of low mass star formation because the final 
outcomes are dramatically different. Low mass star for- 
mation produces isolated or loosely clustered stars while 
high mass star formation yields dense clusters of thou- 
sands of low mass stars alongside relatively few massive 
stars in an ionized nebula. Understanding how high mass 
star formation is initiated and proceeds requires observa- 
tions of the earliest stages of the process, before the newly 
formed massive star has completely disrupted its surround- 
ings. Objects that appear to be high-mass analogues of the 
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Class |Andre. Ward- Thompson fc Barsonvlll99ot ) sources 
observed in low mass s tar formation are now being un- 
cover ed (|Wu et al.ll2005l ; iGarav et al.ll20oi ; iBirkmann et all 
I2007T ). These objects can be studied to determine how close 
an analogue of Class sources they are, enabling the origins 
of the differences between low and high mass star formation. 

G333. 125-0. 562 was orig i nally noted as a potentially 
young object by IGarav et al.l (|2004h and measurem e nts o f 
the mass and temperature of the source bv lLo et alj (|2007l ) 
confirmed the source to be massive dense and cold. It has a 
gas mass of 1.8 x 10 3 Mq, temperature of 13 to 19 K and 
is und etected at wavelengths shorter than 24 /im. iLo et all 
(I2007T ) found G333. 125-0.562 to be a source of SiO emission, 
and the broader width of the SiO line compared to other 
species indicates that this shocked gas tracer species is likely 
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to have been excited in a molecular outflow. This finding is 
interesting in that it suggests that the central source must 
be evolved enough for an organised infall and accretion via 
an envelope (and/or a disk) to have already developed, yet 
before the source is readily detectable in the infrared. 

To further investigate G333. 125-0. 562 we have obtained 
new single dish (CO J = 4 - 3) and interferometric molec- 
ular line (HCO+ J = 1 - 0, SiO J = 1 - and J = 2 - 1) 
data, which are presented in Sj2]and ij3] The integrated emis- 
sion map combined with Spitzer MIPS data reveal that 
G333. 125-0.562 actually contains two components. One of 
these sources dominates the emission and drives the out- 
flow while the nature of the other source is less clear, being 
either an extension of G333. 125-0.562 or perhaps an even 
younger cold core. Maps of molecular line profiles from spa- 
tial positions across the source are presented and these data 
are modelled using a 3D molecular line radiative transfer 
code - MOLLIE (MOLecular Line Explorer; SQJ. The mod- 
elling indicates that heavy depletion of CO and its isotopes 
is taking place, and the abundan ces required for all species 
are consistent with the results of lLo et al.l |2007i ). In >JS]the 
source is compared and contrasted with other sites of high 
mass star formation within the G333 complex. Conclusions 
are drawn in JS] 



calibrator used was Uranus and the secondary calibrator 
was 1646-50 for both 3 and 7 mm observations. The flux 
of Uranus were 7.5 and 1.9 Jy at 3 and 7 mm wavelengths 
respectively. The flux of 1646-50 were 1.6 and 1.9 Jy at 3 and 
7 mm wavelengths respectively. The spectral resolution is 0.8 
km s -1 per channel for both wavelengths. ATCA data were 
reduced with the radio interferometry data reduction pack- 
age miriae0. All ATCA images presented in this work are 
corrected for the primary beam area, the Jansky to Kelvin 
conversion factors are 7.2 and 46.5 K Jy -1 for 3 and 7 mm 
data respectively. 

CO (4 - 3) spectra were obtained with the NANTEN2 
Telescope at Pampa la Bola, Chile. The telescope is 4 me- 
tres in diameter, with angular resolution of 45 arcseconds, 
velocity resolution of 0.6 km s -1 channel - and beam effi- 
ciency of 0.5 at 460-GHz. Observations were made with the 
KOSMA SMART receiver, a dual-frequency 2x8 pixel array 
receiver and the KOSMA array Acousto Optical Spectrom- 
eters (AOS) as backends. 

The 24 and 70 /im images presented in this work are 
part of the Spitzer legacy project 'A 24 and 70 Micron 
Survey of the Inn er Galactic Disk with MIPS' (MIPSGAL; 
ICarev et alJl200Sh . The resolutions are 6 and 18 arcseconds 
at 24 and 70 fxm respectively. The 5cr rms point source 
sensitivities are 1.3 and 73 mjy at 24 and 70 /im. 



2 OBSERVATIONS 

Data presented in this work were collected throughout 2004 
to 2008, comprised of data from the Mopra Radio Telescope, 
the Australia Telescope Compact Array (ATCA)Q and the 
NANTEN2 Observatorjfl Observed molecular transitions 
and continuum emission are listed in Table [T] 

The mapping data were collected between 2004 July 
and 2006 October with the Mopra Telescope, as part of 
the multi- molecular line mapp i ng of the G333 complex (for 
detai ls see iBains et al.1 120061 ; IWong et all 120081 ; lLo et all 
l2009h . The Mopra Telescope is a 22-metre telescope lo- 
cated near Coonabarabran, Australia. It is a centimetre- 
and millimetre-wavelength antenna having a full width to 
half-maximum (FWHM) beam size of ~36 arcseconds and 
velocity resolution of ~0.1 km s -1 channel - at 100-GHz 
|Ladd et al.ll2005l ). The observations were carried out with 
the University of New South Wales Mopra Spectrometer 
(UNSW-MOPS) digital filterbank backend and Monolithic 
Microwave Integrated Circuit (MMIC) receiver, except for 
13 CO and C ls O which were observed with the Mopra Cor- 
relator (MPCOR) and the p revious SIS receive r . For full ob- 
servat i onal setup detail s, see lBains et al. l |2006l ); IWo"n7 et al.l 
|2008l ); lLo et~afl |2009l ). 

High resolution SiO, HCO + and continuum (3 and 7 
mm) imaging of the core were obtained with the ATCA, 
which is an array of six 22-metre antennas at the Paul Wild 
Observatory, Australia. The data was collected during 2007 
October and 2008 March, with array configuration of H75 
and H168, giving a synthesis beam size of 6" x 4" and 5" x 3" 
at 3 and 7 mm wavelengths respectively. For the 3-mm ob- 
servations only five antennas were in use. The primary flux 

1 http://www.atnf.csiro.au 
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3 RESULTS 

3.1 Continuum 

Figure [T] is a map of the ATCA interferometric image of 
3 and 7 mm continuum (thin and thick contours respec- 
tively), overlain on the MIPSGAL 24-/im (a) and 70-/im (b) 
image of G333. 125-0.562. The millimetre continuum emis- 
sion coincides with the 24 and 70 /im emission, while the 
70- /im image reveals that there is another stronger source 
at approximately 5 arcseconds offset to the southwest, with 
fainter 24-/im emission. Herea fter we refer t o the stronger 24- 
/im source as the 'core '. From lJones et al.1 | |200S| ). the core is 
not detectable at 21, 13, 6 and 3 cm wavelengths. Figure[5]is 
a plot of continuum flux-frequency of the core with the 12- 
mm and 3-cm data points being the upper limit. The spec- 
tral index a (S v oc v a , the slope of flux-fre quency plot) cal- 
culat ed using the integrated fluxes at 1.2 (|Mookeriea et all 
|2004 ). 3 and 7 mm is 4.20 ± 0.07. Since these points lie on 
the Rayleigh- Jeans approximation of the Planck function, 
the index should give a good approximation of the actual 
spectral index. The derived spectral index a yields a dust 
opacity spectral slope of /3 = a — 2 = 2.20 ± 0.07. There is 
no clear indication of a turnover frequency (where the opac- 
ity changes from thick to thin) from the flux-frequency plot. 
The implication of this is discussed in iJS] 

3.2 Molecular lines 

Shown in Figure [3] is an ATCA map of the zeroth moment 
(integrated intensity) of the HCO + (1-0) data (thin con- 
tours) and SiO (1-0) data (thick contours), overlain on a 

3 http: / /www. atnf.csiro.au/computing/software/miriad 
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Table 1. List of observed molecula r transitions and continuum of G 333. 125-0, 562 presented in this paper. The columns are: (1) molecule; 
(2) transition; (3) res t frequency i Lovas, Johnson fc Snvd"erlll979h ; (4) upper energy levels from the Cologne Database for Molecular 
Spectroscopy (CDMS. lMiiller et alj200ll . l2005h : (5) calculated critical density: n cr it = A u i/ (a(v)v) where A u i is the Einstein A coefficient 
obtained from CDMS, (c(v)) = 10 — 15 cm 2 is the collision cross section and v is the velocity of collision particles assumed to be 1 km 
s -1 ; (6) observation season; (7) telescope used and (8) the rms noise level. 
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Figure 1. (a) ATCA 3 and 7 mm continuum (thin and thick 
contours respectively) overlain on Spitzer MIPS 24-fim continuum 
image of G333. 125—0.562. The contour levels are at 0.03 and 0.06 
Jy beam -1 for the 3-mm continuum, 0.002 and 0.003 Jy beam -1 
for the 7-mm continuum, (b) Same contoured data but overlaid 
on MIPS 70-/im image. The crossed ellipse at lower left corner 
represent the beam size of the 3 and 7 continuum data. 




frequency [GHz] 



Figure 2. A plot of continu um peak flux-frequency of the core. 
The 3-cm lljones et al.ll2008r i and 12-mm contin uum data points 
are th e upper limit, the 1.2-mm data is from Mookcrjca ct al. 
(2001). The plot shows a spectral index a of 4.2 for the three 
millimetre data points (dashed line) and no clear sign of turning 
point. 



MIPS 24-/im image (grey scale) of the core. The elongated 
shape of the HCO + emission comes from the blue-shifted line 
wings of molecular outflow as shown in Figure [4^,, an inte- 
grated emission map over a velocity range of —64 to —70 km 
s -1 . The ATCA HCO + spectra at the peak position (Figure 
|4p) clearly shows the blue-shifted line wing not seen in the 
Mopra HCO + spectra (Figure 2Ji). No red component of the 
outflow is seen in the ATCA data, which is likely spatially 
oriented such that it is outside the primary beam. Further 
details are presented in Section[3]on the modelling of HCO + 
lines. 

Several 24-/xm sources are clearly present yet the Spitzer 
70- /im emission (Figure [T]) is confined to the brighter source 
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Figure 3. Integrated intensity map of HCO" 1 " (1 - 0; thin con- 
tours) and SiO (2 — 1; thick contours) overlain on Spitzer MIPS 
24-/im image (grey scale) of G333. 125-0. 562. Integrated velocity 
range is from -80 to -50 km s" 1 for both HCO+ and SiO. The 
lowest contour level is 25 K km s _1 with increments of 5 K km 
s _1 for HCO + and 20 K km s^ 1 with increments of 5 K km s _1 
for SiO. The cross marks the position of the 3 and 7 mm con- 
tinuum peaks, and the crossed ellipse represents the beam size of 
the molecular line data. The two distinctive molecular emission 
peaks are denoted by 'A' and 'B'. 



to the northwest, marked 'A' in Figure [3] The southeast 
HCO + emission is marked 'B'. Furthermore, the SiO (2 - 1) 
integrated intensity data (thick contours) is also confined to 
the northwest sourc e. Based on its a ssociation with the SiO 
emission detected bv lLo et all (|2007f l and the spectral energy 
distribution, the northwest source ('A') was identified as the 
primary source in this study. The nature of the secondary 
source and its possible association with G333. 125-0.562 is 
discussed in Section [5] 

Note that the HCO + spectrum does not show clear signs 
of infall (asymmetry profile) unlike the lower resolution Mo- 
pra data; the implication of this is discussed in Section [5] 
Figure [5] [6] and [7] show the channel maps of the ATCA 
HCO + and SiO data respectively. 

Figure [8] show the line profile of CO (4 - 3) obtained 
with the NANTEN2 telescope, with a beam size of 45 arc- 
seconds. The line profile shows two emission peaks at —59 
and —50 km s _1 where the G333 cloud itself contributes to 
the —50 km s" 1 peak. Blue-shifted emission is also present, 
from ~ —65 to ~ —80 km s _1 , consistent with the blue 
shifted emission in the ATCA spectrum (Figure [4]), while 
the red-shifted component is confused with the —50 km s _1 
ambient cloud emission. CO (4 - 3) is a moderate excitation 
line that should better trace the warmer central regions of 
G333. 125-0.562 than lower excitation lines such as CO (1 
- 0) which typically contain contributions from low density 
emission along the line of sight through the molecular cloud 
complex. 

Finally, using new and existing Mopra and ATCA data, 
spatial grids of HCO + line profiles from positions across the 
face of G333. 125-0.562 were produced. These line profiles 
encode information about the velocity, density, turbulence 
and temperature fields inside the core. In the next section, 
a fully 3D molecular line profile radiative transfer code is 
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Figure 4. (a) Integrated emission map of the ATCA measure- 
ment of the HCO+ line towards source 'A', over a velocity range 
of —64 to —70 km s~ 1 , showing the blue-shifted line wing con- 
tributing to the outflows that is not seen in the Mopra data. The 
contour levels are at 3.3, 6.6 and 9.9 K beam . (b) Line pro- 
file at the peak HCO+ (ATCA) position showing the blue-shifted 
line wing, (c) HCO + spectra obtained with Mopra showing the 
absence of blue-shifted line wing. The dashed lines highlight the 
blue-shifted line wings velocity range. 

used to closely model this set of line profiles and so the data 
is presented alongside the model line profiles. 



4 RADIATIVE TRANSFER MODELLING 

To verify and extend the observational analysis, a 3D molec- 
ular line transport radiative transfer code - MOLLIE was 
used to model the primary source in G333. 125-0.562. MOL- 
LIE takes in a 3D grid of physical parameters (e.g. den- 
sity, temperature, velocity, turbulence, chemical abundance) 
of a source, then simulates the molecular emission pro- 
files for comparison with observed line profiles. Example s 
of the use of the code can be found in iKeto et all (|200-4h ; 
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Figure 5. Channel maps of ATCA HCO" 1 " (1 - 0) of the core. The lowest contour level is at 1.0 K with increments of 1.0 K. The cross 
marks the position of 3 and 7 mm continuum, and the crossed ellipse represents the beam size of the molecular line data. 



Redman et alj (|2006l ) ; ICarolan et al. | (|200Sl , l2009D . A model 
of the source was constructed based on the observational 
parameters (chemical abundance, density, v e locity , gas tem- 
perature and line intensity) from lLo et all (|2007h and ref- 
erences therein, with the aim of reproducing line profiles of 
13 CO (1 - 0), C 18 (1 - 0), CS (2 - 1) and HCO+ (1 - 
0), and using the NANTEN2 CO (4 - 3) transition as a 
constraint to the model. 

We have generated over a thousand models comprised 
of different physical parameters that give rise to promi- 
nent features in the line profiles: density and abundance 
(line strength), temperature (line strength and line widths), 
infall/expansion velocity (asymmetric self-absorption line), 
outflow velocity (line wings). By fitting multiple molecular 



transition line profiles, the number of possible best fit mod- 
els is greatly constrained, leaving only a handful to choose 
from. The physical parameters of these limited number of 
possible models vary within a small threshold (a factor of 
two). The modelling consists of two parts. The first reduces 
the number of possible models by fitting the line profiles of 
multiple Mopra molecular lines. In the second part we refine 
the model with high resolution ATCA data. Details of the 
modelling is described in the following paragraphs. 

The radius of the core was taken to be R = 0.5 pc. The 
temperature increases linearly from the edge of the core (12 
K) to the centre (17 K). The density distribution is modelled 
as a simple Plummer sphere p = (p c Rp) / (R^ + r 2 ) such that 
the density is approximately constant in the centre of the 
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Figure 6. Channel maps of ATCA SiO (1 - 0) of the core. Contour level starts at 1.0 K with increments of 1.0 K. The cross marks the 
position of 3 and 7 mm continuum, and the crossed ellipse represents the beam size of the molecular line data. 



core at p c — 5 x 10 5 cm~ 3 and then falls off as r -2 . The 
scaling factor of R p = 0.1R is quite small which means that 
the density distribution is close to a r~ power law. This 
density distribution was chosen to be consistent with the 
fall off in the strength of the line profiles away from the 
centre of the core, as well as with the dust emission profile. 



4.1 Modelling of Mopra data 

Typical line widths for species other than SiO are 4-5 km 
s _1 . Thermal broadening (due to gas at 12 - 17 K) is far 
too small to generate these line widths which implies that 
dynamical processes such as turbulence, infall, outflow and 
rotation must be responsible for line widths. Accordingly, 



the model also incorporates a directed molecular outflow, 
radial infall and bulk turbulence. Infall is required because 
lines such as HCO + (Mopra data) exhibit a self-absorption 
profile with a stronger blue-shifted peak to that observed 
in nearby low mass infall sources. The line split due to in- 
fall is marginally seen in the asymmetrical shape of the line 
profile, in which the blue shifted emission is stronger due 
to self- absorption of the red emission (Figure O, suggesting 
the infall velocity is comparable to the turbulent velocity. 
The self-absorption of the red shifted emission in this way, 
rather than widely split or blended profiles, is routinely seen 
in single dish observations of infall in low mass star form- 
ing clouds, and is seen here with the larger line widths and 
splitting of a massive star forming region. Isotopes such as 
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Figure 7. Channel maps of ATCA SiO (2 - 1) of the core. Contour level starts at 1.5 K with increments of 1.5 K. The cross marks the 
position of 3 and 7 mm continuum, and the crossed ellipse represents the beam size of the molecular line data. 



H 13 CO + show single peaked Gaussian shape profiles and 
turbulence is required to give a structureless broad line pro- 
file at these temperatures. Turbulent velocities of 1.5 km 
s _1 and infall velocities of —1.5 km s _1 provided the best 
fit to the lines. The addition of a low speed molecular out- 
flow (velocity a few times the infall speed) is included for 
completeness but is in fact not constrained well by the ob- 
servational data: the Mopra beam is comparable to the size 
of the entire core and the outflow material will thus only 
occupy a small fraction of the beam, as evident from the 
presence of outflow line wings in the ATCA HCO + data 
(Figured]). 

The chemical abundances used were tho s e deri ved from 
the analysis of the observations in lLo et all (|2007h . For all 



molecules depletion, due to the freeze-out the molecules onto 
dust grains, is required in order to fit the data. However, 
the degree of depletion required varies widely between the 
species. 13 CO and C ls O require a reduction of an order of 
magnitude of their gas phase abundance within r < 0.9R, 
from 1 x 1(T 5 and 8 x 1(T 7 to 1 x 1CT 6 and 8 x 1CT 8 for 
13 C0 and C 18 respectively (Table©. This is entirely con- 
sistent with the presence of N 2 H + in the centre of the core 
which should be destroyed by CO and is also consistent with 
integrated intensity maps of these species which show the 
CO isotopes are only weakly present above the background 
molecular cloud. More modest depletion is implied for the 
HCO + data and the modelling suggests that this species is 
depleted within r < 0.3R. Little freeze-out of CS is implied 



8 Lo et al. 




-60 
Velocity (km/s) 



Figure 8. CO (4 - 3) line profile of the core obtained with NAN- 
TEN2. The line profile shows two emission peaks at —59 and —50 
km s _1 . 



by the modelling which calls for a depletion only within the 
central regions at r < 0.1R. The core (envelope) and outflow 
parameters of the best fit model are summarised in Table [2] 
the modelled spectra are shown in Figure [9] 



4.2 Modelling of the ATCA data 

Based on the best fit model of Mopra data, we further re- 
fined the input parameters to replicate the blue-shifted line 
wings found in the ATCA HCO + line profiles. We adopted 
an 'hourglass'-shaped outflow model (|Rawlings et alj|2004f ) 
with an opening angle of 50 degrees. In order to reproduce 
the centre line intensity, the core density has to be doubled 
from the Mopra data model (Table [2}. The outflow density 
is also an order of magnitude higher and temperature of 50 
K to give the line wings intensity, this is consistent with the 
suggestion that the Mopra HCO + data suffers from beam 
dilution as the source does not fill up the beam, and that 
the interferometric HCO + data has filtered out the extended 
structure of the core. The extent of the HCO + line wings are 
reproduced with an outflow velocity of 12 km s _1 , while a 
turbulent velocity of 3.0 km s _1 gives the smooth blending 
profile between the centre line component and the outflow 
line wings; the outflow is oriented such that the blue-shifted 
outflow lobe is south-west of the source and the red-shifted 
lobe falls outside of the primary beam. The HCO + abun- 
dance has to be increased by two orders of magnitude in 
order to reproduce the observed intens i ty. T his is consis- 
tent with the findings of iRawlings et al.l \2Q04 i that HCO+ 
abundance is enhanced in outflows. Shown in Figure [Till are 
sample spectra of the observed ATCA HCO + overlain on 
the model spectra. 

In summary, a self-consistent radiative transfer model 
was generated using t he size, t e mpera ture, density and 
abundances reported in iLo et all (|2007r i. Supersonic infall 
and turbulence and heavy depletion of CO, due to freeze- 
out, is required to produce a good fit to the line widths, 
shapes and strengths. 



5 DISCUSSION 

G333. 125-0. 562 is located in the G333 giant molecular cloud 
complex which has been the subject of a mult i- molecular line 
survey project carried out by the authors. Several papers of 
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Velocity (km/s) 

Figure 9. Modelled line profiles (thin lines) and Mopra observa- 
tions line profiles (thick lines) of the core. We aim to reproduce 
the line strength and line width of the observed molecular tran- 
sitions. Note the observed emissions at ~ —50 km s -1 is due 
to ambient cloud emissions and thus not present in the model 
spectra. 



resul t s have now been published from t his work (|f3ains et al.l 
120061 ; IWong et al.ll2008l ; ILo et al.ll2009t ) and so it is now pos- 
sible to compare G333. 125-0.562 with other sites of high 
mass star formation within G333. 

Models of high mass star formation emphasise the role 
of turbulence (|McKee fc Ostri kcr 2007), so quantifying the 
variation in the properties of turbulence throughout G333 



Modelling 



of massive dense cold core G333. 125-0. 562 9 



Table 2. Summary of characteristic model parameters for the 
ATCA and Mopra data determined using the MOLLIE radiative 
transfer code. All parameters vary with functional forms that are 
described in detail in the text. Density values listed here are the 
central densities following a density distribution of a simple Plum- 
mer sphere as detailed in Section [4] The listed values are uncer- 
tain by a factor of two, as derived from the number of plausible 
models. 



Mopra HCO+, CS, C ls O and "CO 


Envelope 


Outflow 




+ Core^ 




Density (xlO 5 cm -3 ) 


5.0 


0.1 


Velocity (km s _1 ) 


-1.5* 


3.0 


Turbulent width (km s — 1 ) 


1.5 


1.5 


Temperature (K) 


17 


17 


Abundance C ls O/H 2 (xlCT 7 ) 


8.0 


8.0 


Abundance 13 CO/H 2 (xlCT 5 ) 


1.0 


1.0 


Abundance CS/H 2 (xlCT 10 ) 


1.0 


1.0 


Abundance HCO+/H 2 (xlCT 10 ) 


6.8 


6.8 


ATCA HCO+ 


Core 


Outflow 


Density (xlO 5 cm -3 ) 


10.0 


3.0 


Velocity (km s _1 ) 


-1.5" 


12.0 


Turbulent width (km s" 1 ) 


3.0 


3.0 


Temperature (K) 


17 


50 


Abundance HCO+/H 2 (xlCT 10 ) 


6.8 


34.0 



a Mopra beam resolution covers both the core and envelope. 
* Infall velocity. 
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Figure 10. Modelled HCO + spectra (thin lines) overlaid on the 
observed spectra (ATCA; thick lines) across the core. The num- 
bers in top left brackets are offsets from the HCO + peak in arc- 
seconds. 



is important to the testing of such models. Turbulent veloc- 
ities in G333. 125-0.562 are highly supersonic and large by 
low mass star formation standards. However, they are signif- 
icantly less than those measured in the three most massive 
star forming regions in G333 (~ 6 km s" 1 ; Bains et al, 
in prep.). The core has a lower mass and a lower turbulent 
(and infall) velocity that might be expected simply from 



virial theorem ar guments, where th e mass of the cold core 
is 2.2 x 10 3 Mq (|Garav et alj|2004 ). an order of magnitude 
lower than the three most massive sources in G333. 

The coincidence of the continuum sources at four dif- 
ferent wavelengths (3 and 7 mm, 24 and 70 /jm) further 
confirms that the internal heating source is deeply embed- 
ded in the cold core. Whether the continuum emission is due 
to warm dust or free-free emission from a deeply embedded 
Hn region can be examined using the spectral index a (see 
Section [3]). The three millimetre wavelength (1.2, 3 and 7 
mm) data points yield a spectral index a of 4.2 ± 0.07 and 
dust emissivity /3 of 2.2, w hich is comparab le to the dust 
emissivity of j3 — 2 found bv lHill et all (|2006l ) for cold cores. 
Thus, to at least 7-mm wavelength, the continuum emission 
is due to thermal dust emission rather than free-free emis- 
sion. At longer wavelengths, there is no detectable (free- free) 
continuum emission from the core, thus it is possible that 
an Hn region still has not yet developed at this stage. 

The secondary molecular source (denoted 'B' in Figure 
[3]) has no detectable millimetre, 24 nor 70 /im continuum 
emission, thus no indication of an internal heating source. 
We speculate this source could be an even less evolved ob- 
ject. 

Contrary to the Mopra HCO + spectra profiles, the high 
resolution ATCA data does not show infall asymmetry. This 
is likely due to extended emission, such as the contribution 
from the envelope, where infall takes place, is filtered out in 
the interferometry data. 

Two possible explanations for the absence of outflow 
line wings in Mopra HCO + line profiles are lack of sensitivity 
and/or the outflow (spatial) scale. The ATCA HCO+ data 
has a spatial resolution of 6 x 3 square arcseconds while the 
Mopra data has resolution of 36 arcseconds, if the outflow is 
not widely spread and only partially fills the Mopra beam, 
then it is heavily diluted. Thus line wings are not seen in 
the Mopra spectra profiles. An abundance enhancement is 
also supported by the modelling, which requires a factor 
of 5 larger HCO + abundance for the outflow component to 
reproduce well the higher resolution ATCA data than it does 
for the Mopra data. 



6 CONCLUSIONS 

G333. 125-0.562 is an excellent example of the very early 
stages of massive star formation. The cloud is cold, moder- 
ately turbulent and possesses an outflow which indicates the 
presence of a central driving source. HCO + is enhanced in 
the outflow. The spectral index of the primary source indi- 
cates the continuum emission at wavelengths shorter than 
7-mm is due to thermal dust emission. Undetectable contin- 
uum emissions at longer wavelengths suggests the absence 
of free- free emission from ionised gas, thus an Hn region has 
not yet developed. We have modelled the molecular emission 
using the 3D radiative transfer code MOLLIE, obtaining the 
best match to the data with a model where there is a mas- 
sive, dense core containing an embedded source driving an 
outflow, but without the presence of detectable ionized gas 
from an Hn region. 



10 Lo et al. 



ACKNOWLEDGMENTS 

NL acknowledges partial support from Center of Excel- 
lence in Astrophysics and Associated Technologies (PFB 
06) and Centra de Astrofi'sica FONDAP 15010003. NL's 
postdoctoral position at CEA/Irfu was funded by the Ile- 
de-France Region. MPR acknowledges support from a Sci- 
ence Foundation Ireland Research Frontiers grant. PAJ ac- 
knowledges partial support from Centro de Astrofi'sica FON- 
DAP 15010003 and the GEMINI-CONICYT FUND. This 
work was supported by a Science Foundation Ireland Re- 
search Frontiers grant. The Mopra Telescope and ATCA 
are part of the Australia Telescope and are funded by the 
Commonwealth of Australia for operation as National Fa- 
cility managed by CSIRO. The UNSW-MOPS Digital Filter 
Bank used for the observations with the Mopra Telescope 
was provided with support from the University of New South 
Wales, Monash University, University of Sydney and Aus- 
tralian Research Council. The NANTEN2 Observatory is a 
collaboration between Nagoya University, Osaka University, 
KOSMA, Universitat zu Koln, Argelander-Institet Univer- 
sitat Bonn, Seoul National University, ETH Zurich, Univer- 
sity of New South Wales and Universidad de Chile. This 
research has made use of the NASA/IPAC Infrared Science 
Archive which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology, under contract with 
NASA. 



Lovas F. J., Johnson D. R., Snyder L. E., 1979, ApJ Supp, 
41, 451 

McKee C. F., Ostriker E. O, 2007, ARA&A, 45, 565 
Mookerjea B., Kramer O, Nielbock M., Nyman L.-A., 

2004, A&A, 426, 119 

Midler H. S. P., Schloder F., Stutzki J., Winnewisser O, 

2005, J. Molecular Structure, 742, 215 

Midler H. S. P., Thorwirth S., Roth D. A., Winnewisser 

G., 2001, A&A, 370, L49 
Rawlings J. M. O, Redman M. P., Keto E., Williams D. A., 

2004, MNRAS, 351, 1054 

Redman M. P., Keto E., Rawlings J. M. O, 2006, MNRAS, 
370, LI 

Wong T., Ladd E. F., Brisbin D., Burton M. O, Bains 
I., Cunningham M. R., Lo N., Jones P. A., et al., 2008, 
MNRAS, 386, 1069 

Wu Y., Zhu M., Wei Y., Xu D., Zhang Q., Fiege J. D., 

2005, ApJ Lett, 628, L57 

This paper has been typeset from a TgX/ F/IgX file prepared 
by the author. 



REFERENCES 

Andre P., Ward-Thompson D., Barsony M., 1993, ApJ, 
406, 122 

Bains I., Wong T., Cunningham M., Sparks P., Brisbin D., 
Calisse P., Dempsey J. T., Deragopian C, et al., 2006, 
MNRAS, 367, 1609 

Birkmann S. M., Krause O., Hennemann M., Henning T., 
Steinacker J., Lemke D., 2007, A&A, 474, 883 

Carey S. J., Noriega- Crespo A., Mizuno D. R., Shenoy S., 
Paladini R., Kraemer K. E., Price S. D., Flagey N., et al., 
2009, PASP, 121, 76 

Carolan P. B., Khanzadyan T., Redman M. P., Thompson 
M. A., Jones P. A., Cunningham M. R., Loughnane R. M., 
Bains I., et al., 2009, MNRAS, pp 1353 — h 

Carolan P. B., Redman M. P., Keto E., Rawlings J. M. O, 
2008, MNRAS, 383, 705 

Garay G., Faundez S., Mardones D., Bronfman L., Chini 
R., Nyman L.-A., 2004, ApJ, 610, 313 

Hill T., Thompson M. A., Burton M. G., Walsh A. J., 
Minier V., Cunningham M. R., Pierce-Price D., 2006, MN- 
RAS, 368, 1223 

Jones D. I., Protheroe R. J., Crocker R. M., 2008, Publi- 
cations of the Astronomical Society of Australia, 25, 161 

Keto E., Rybicki G. B., Bergin E. A., Plume R., 2004, ApJ, 
613, 355 

Ladd N., Purcell C, Wong T., Robertson S., 2005, Publi- 
cations of the Astronomical Society of Australia, 22, 62 

Lo N., Cunningham M., Bains I., Burton M. O, Garay G., 
2007, MNRAS, 381, L30 

Lo N., Cunningham M. R., Jones P. A., Bains I., Burton 
M. G., Wong T., Muller E., Kramer C, et al., 2009, MN- 
RAS, 395, 1021 



